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The global variables like the transverse energy at midrapidity, the charged particle multiplicity 
at midrapidity and the total multiplicity of charged particles are evaluated in the single-freeze-out 
\Q ■ statistical model for different centrality bins at RHIC at ^/snn = 130 and 200 GeV. Full description 

f"^ ' of decays of hadron resonances is applied in these estimations. The geometric parameters of the 

model are obtained from the fit to the final data on the pr spectra. The predicted values of the 
£SJ . global variables agree qualitatively well with the experimental data. The centrality independence 

, | of the total number of charged particles per participant pair has been also reproduced. 
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PACS numbers: 25.75.-q, 25.75.Dw, 24.10.Pa, 24.10.Jv 



(N) ■ I. INTRODUCTION 

> . 

In the previous paper pj the extensive analysis of the centrality dependence of two measured global variables, 
transverse energy (cIEt / dr)\ m id) and charged particle multiplicity {dN c h / dw\ m id) densities at mid-rapidity as well as 
their ratio, was delivered. The analysis was done in the framework of the single freeze-out statistical model 
for the Au-Au collisions at RHIC at ^Jsnn = 130 and 200 GeV. The main idea of this method is as follows. Thermal 
and geometric parameters of the model are established from fits to particle yield ratios and pt spectra, respectively. 
Then, with the use of these parameters both densities dET/dr) and dN c h/dn as well as the total charged-particle 
multiplicity can be estimated and compared with the data. The main reason for doing it is that the transverse energy 
and charged particle multiplicity measurements are independent of hadron spectroscopy (in particular, no particle 
identification is necessary), therefore they could be used as an additional test of the self-consistency of a statistical 
model. In the present paper the above-mentioned program will be performed once more for the reasons stated below. 

First of all fits of the geometric parameters [5j were done with the use of the preliminary data for the pr spectra 
measured at RHIC at ^snn = 200 GeV 0,0- But the final data 0,0] differ substantially from the preliminary ones. 
Additionally the data points for the \ 2 analysis were digitized from the plots. Also not all centrality bins were fitted 
in 0, . All of these resulted in the rough qualitative agreement of the model predictions for dN c h / drj and dEr / drj 
for RHIC at ^/snn = 200 GeV 1.]. And the centrality independence of the total multiplicity of charged particles 
per participant pair was not reproduced. The evidence for such scaling of the total multiplicity was reported by the 
PHOBOS Collaboration Therefore it was very tempting to check whether the \ 2 analysis of the final data for 
the pt spectra of identified charged hadrons would improve the predictions for transverse energy and charged particle 
multiplicity densities and the total charged-particle multiplicity. Such a check is nontrivial and does not mean simply 
adding the integrals of the experimental momentum distributions because: (a) spectra are measured in the limited 
ranges of px and the very important low-px part is not detected, (b) spectra of stable neutral particles like neutrons 
and K® L are not measured generally. 

As a statistical model the single freeze-out model is taken (for details see Fiol]). The model succeeded in the accurate 
description of ratios and pr spectra of particles measured at RHIC 0, H, Q|- The main postulate of the model is the 
simultaneous occurrence of chemical and thermal freeze-outs, which means that the possible elastic interactions after 
the chemical freeze-out are neglected. The conditions for the freeze-out are expressed by values of two independent 
thermal parameters: T and fig. The strangeness chemical potential (is is determined from the requirement that the 
overall strangeness equals zero. The second basic feature of the model is the complete treatment of resonance decays. 
This means that the final distribution of a given particle consists not only of the thermal part but also of contributions 
from all possible decays and cascades. Feeding from week decays is included as well. 
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II. THE FOUNDATIONS OF THE SINGLE-FREEZE-OUT MODEL 

The main assumptions of the model are as follows: (a) the chemical and thermal freeze-outs take place simultane- 
ously, (6) all confirmed resonances up to a mass of 2 GeV from the Particle Data Tables |l2| are taken into account, 
(c) a freeze-out hypersurface is defined by the equation 



r ' 2 — const , (1) 
(d) the four- velocity of an element of the freeze-out hypersurface is proportional to its coordinate 



«" = - = - (l,T'T'T) ' ( 2 ) 

T TV t t t J 

(e) the following parameterization of the hypersurface is chosen: 

t = rcoshaii coshaj_, r x = r sinh a_i_ cos <^>, r B = r sinh ax sin r 2 = rsinhan coshaj_, (3) 
where a\\ is the rapidity of the element, oc\\ = tanh - (r z /t), and a± controls the transverse radius: 



+ r 2 = r sinh a± < p max , (4) 

where the restriction on the transverse size has been introduced, so p ma x gives the maximal transverse extension 
of the gas in the central slice during the freeze-out. This means that two new parameters of the model have been 
introduced, i.e. r and p max , which are connected with the geometry of the freeze-out hypersurface. In general a 
freeze-out addresses a 6-dimensional phase-space but the assumption of the Hubble-like expansion (the velocity at 
the freeze-out is proportional to the coordinate), Eq. J2J), reduces the number of independent variables to 3. And 
the assumption of the cylindrical symmetry further reduces this number to 2. Another point is how realistic the 
assumption of the Hubble flow is. Recent results of Refs. [l3l[T3 | indicate that the Hubble flow of the form of Eq. (J2J 
can develop in 130 and 200 GeV Au-Au collisions within times 6-15 fm. This is roughly the scale of the freeze-out 
initializing time predicted in this analysis for all but peripheral bins (see Table QJ. 

From Eq. one can see that the beginning of the freeze-out process starts at t^P Q = r and r = in the c.m.s., 
which is also the laboratory frame in the RHIC case. At this moment the volume of the gas can be estimated as 

V#> = 2irrp 2 max , (5) 

which is simply the volume of a tube with a length 2r and a radius p max {%t is the maximal possible extension of 
the gas in the longitudinal direction at tj Q ). In the central slice the freeze-out ceases at vf' a = \Jt 2 4- p^ a;E and it 
takes place at p = p max . This does not mean the end of the freeze-out process completely. It is the moment when 
the gas decouples into two pieces symmetrical with respect to the plane 2 = and from then the freeze-out proceeds 
into two opposite directions of the collision axis z. For the considered hypersurface, Eq. ft}, this process extends to 
infinity but to be more realistic a limiting condition for ay will be assumed in Sec. IIV CI namely a maximal possible 
value of a || called a™ ax will be postulated. Then the moment when the freeze-out ceases at all will be equal to 

tfl=co S hal a *-tfl. (6) 

(2) 

At t j the volume of the region where the freeze-out process has just taken place or is to happen can be estimated as 

Vf 2 l = 27Tv/t 2 + p 2 max p 2 max . (7) 
The transverse velocity in the central slice can be expressed as a function of the transverse radius 

p ± (p) = tanh a ± = P . (8) 
i/r + p l 
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The maximum value of /3± called the maximum transverse-flow parameter (or the surface velocity) is given by 

pmax _ Pmax _ Prnax / T ^ 

V T ' 2 + Pmax V 1 + (Prnax/T) 2 

The invariant distribution of the measured particles of species i has the form 0, Q 



d 2 p T dy 



p"do), /<(?•«), (10) 



where da^ is the normal vector on a freeze-out hypersurface, p ■ u = p^u^ , u M is the four-velocity of a fluid element 
and fi is the final momentum distribution of the particle in question. The final distribution means here that Ji is the 
sum of primordial and simple and sequential decay contributions to the particle distribution (for details see flcLllfijl). 
With the use of Eqs. (J2J and J3J), the invariant distribution ljTU|) takes the following form: 



+ 00 Pmax /t 27T 

dNi 



d 2 pr dy 



r 3 



J da || J sinhaj_d(sinhax) / d£ (p ■ u) fi(p ■ u) , (11) 



-oo 

where 

p ■ u — rriT cosh (an — y) cosha^ — px cos£ sinhaj_ . (12) 

After changing the integration variable an to ay + y the invariant distribution (|ll|l still has the same form but with 
p ■ u replaced now by p^u: 

p ■ u = rax cosh ay cosh a± — pt cos £ sinh aj_ . (13) 
II . If a || is limited and has its maximal value a™ 



(fTU|l reads 



dNi 3 

— T 



dan / sinhaj_d(sinhaj_) / d£ (p ■ u) fi(p ■ u) — I da (p ■ u) fi(p ■ u) . (14) 



o o 



d 2 p T dy 

The same changing of the integration variable an — > an + y can be done and one obtains 



dN 3 

— T 



d 2 pr dy 



J cZct | J sinhaxd(sinhax) / d£ (p ■ u) fi(p ■ u) . (15) 



The distribution l|15fl is still boost invariant because rapidity differences are invariant with respect to the longitudinal 
boosts. 

In practical, i.e. numerical, calculations of the invariant distribution given by Eqs. (|llf) and l|13(l a cut in the integral 
variable a\\ has to be done. So for the numerical reasons the range of a\\ is always finite and — a| ut < a\\ < a^ ut . The 

value of a| u * is chosen in such a way that the remaining "tails" of the integral over a\\ are negligible. It has turned out 
that it is enough to put ajj nt = 3.0 to determine the geometric parameters r and p max exact to four decimal places. 

It should be reminded here that in this paper fits are done to the spectra measured at midrapidity (y = 0). Also 
the estimates for dEx/dr] and dN c h/dr] are done for y = 0. What is more, the values of a™ ax predicted in Sec. II V 01 
are greater than 3.0. Thus for the practical calculations at y = 0, the case with unlimited ay covers the case with the 
finite range of an. 
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III. TRANSVERSE ENERGY AND CHARGED PARTICLE MULTIPLICITY DENSITIES 

The experimentally measured transverse energy is denned as 



suit 



(16) 



where 9i is the polar angle, denotes E t — mjv (tojv means the nucleon mass) for baryons, Ei + for antibaryons 
and the total energy Ei for all other particles, and the sum is taken over all L emitted particles [lfil Il7j . 
The pseudorapidity density of particle species i is given by 



dNj = r d2 dy dNj 
drj J dr] d 2 px dy 



= / d 2 PT 



P 



dNi 



Ei d 2 p T dy ' 

Analogously, the transverse energy pseudorapidity density for the same species can be written as 



(17) 



dE T ,i 



d 2 p T Ei 



p T dy dNi 



d Pt Pt 



Et dNi 



dr) J r± 1 p drj d 2 pT dy J Ei d 2 px dy 

For the quantities at midrapidity one has (in the c.m.s., which is the RHIC case) 



(18) 



dNj 
drj 



mid 



d pt 



Pt 



dN, 



rriT d 2 pT dy ' 



(19) 



dE T ., 
dy 



[\d 2 pTPT^^J^ Ty ^ = baryon 



mid 



= < J d 2 p T Pt 



J d 2 p T Pt 



tut d' 2 pT dy ' ' 



antibaryon 



dNj , 
d 2 pT dy ' 



others . 



The overall charged particle and transverse energy densities can be expressed as 



(20) 



dNch 
dr) 



drj 



(21) 



dE T 
drj 



raid 



. - dE_ 



TA 



mid 



(22) 



where A and B {B C A) denote sets of species of finally detected particles, namely the set of charged particles 
B = {7r + , 7r~, K + , K~, p, p}, whereas A also includes photons, K {) L . n and n [T^. 



IV. RESULTS 



A. Determination of geometric parameters 



Analyses of the particle ratios and pt spectra at various centralities in the framework of the single freeze-out model 
were done for the preliminary RHIC data at ^snn = 200 GeV |6j, |jj in [5|- This approach proceeds in two steps. 
First, thermal parameters T and /j.b are fitted with the use of the experimental ratios of hadron multiplicities at 
midrapidity. It is assumed that these values are independent of the centrality. This is reasonable since the very weak 
centrality dependence of the particle ratios has been observed so far. Recent works of Ref. [n| 1 2(1 l2l| have just 
confirmed this assumption. After then two next parameters, r and p ma x, are determined from the simultaneous fit 
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to the transverse- momentum spectra of n , K j) and p. Both fits are performed with the help of the \ 2 method. 
Since the preliminary data for the pr spectra @, Q differ from the final data 0, || and not all bins were fitted in |{| , 
the fit procedure for determination of the geometric parameters of the model, r and Pmax, has been performed again 
for the purposes of this paper. Additionally, the PHENIX case at ^/snn — 130 GeV has been worked out once more, 
since the first published data were for three bins only |2^ . The data for the next two bins, which were not fitted in 
[To|. were added in the later report ^23jj. Also the BRAHMS spectra at various centralities for the measurement at 
^/snn = 200 GeV [24j have not been fitted within the single-freeze-out model until now. The data points with pt > 3 
GeV have been excluded in the \ 2 analysis. The thermal parameters for the three cases of the maximal RHIC collision 
energy have been taken from the newer studies of the particle abundance ratios [2(1 l2l| . According to these results, 
the thermal parameters are different for two experiments, PHENIX and STAR, in spite of the same collision energy. 
Fits to mostly PHENIX ratios gave the lower freeze-out temperature [2(j than fits to STAR data only [2*l|. This 
observation has been confirmed in the present analysis: fits to the PHENIX spectra with substitution of the STAR 
freeze-out thermal parameters (T = 160.0 MeV and ps = 24.0 MeV [2l|) and fits to the STAR spectra with the use 
of the PHENIX freeze-out thermal parameters (T = 155.2 MeV and ps = 26.4 MeV [2(J) have proven to be much 
worse than the fits done within the same experiment and finally reported in this work (here " worse fit" means that its 
X 2 /NDF is greater). Also it has turned out that BRAHMS spectra fit better with the PHENIX thermal parameters. 
The older estimates of thermal parameters for RHIC at ^/s^n = 200 GeV (T = 165.6 MeV and p B = 28.5 MeV) 
given in Ref. [f| have proven to lead to fits to PHENIX spectra which are not statistically significant for all centralities, 
X 2 /NDF > 1.3. It should be also mentioned at this point, that the values of thermal parameters (T = 177 MeV and 
Pb = 29 MeV) obtained in Ref. |26j result in fits which are much worse, x 2 /NDF s» 3. 

The final results for the geometric parameters p m ax and r are gathered in Tabic [I] together with the corresponding 
values of x 2 /NDF for each centrality class additionally characterized by the number of participants N part . Note that 
besides the most peripheral bins of the PHENIX measurements all fits are statistically significant. Other physical 
quantities like the surface velocity f3™ ax , the volume at the beginning of the freeze-out K 1 ], the maximal freeze-out 

time at the central slice and the corresponding volume vj 2 ] , which have been defined in Sec. [HJ are also given 
there. Note that values of p max and r fitted for the same spectra as in (the — 5%, 15 — 30% and 60 — 92% 
centrality bins for PHENIX at ^snn = 130 GeV) are the same with 5% accuracy as given in JlJ}, besides the 
peripheral bin. The model calculations take into account full feeding from weak decays besides the PHENIX case at 
ytsNN = 200 GeV where protons (antiprotons) from A (A) decays are excluded. This needs some comments because 
the STAR Collaboration claims that its pion spectra are corrected for weak decays ,9]. However, in the present 
analysis the best quality of the fit expressed by the value of % 2 /NDF in the STAR case is for full feeding from weak 
decays. Namely, if pions from A (A) decays are subtracted then x 2 /NDF = 0.37 for the — 5% centrality bin, whereas 
if pions from all possible weak decays are excluded % 2 /NDF = 0.88 for the same bin. With the full inclusion of pions 
from weak decays % 2 /NDF = 0.29 (see Table|IJ. That the single-freeze-out model with all the resonances and the full 
feeding from weak decays can fit the STAR spectra has been already noticed in |9j but, as it is stated there, with not 
satisfactory % 2 /NDF. 

Values of the geometric parameters p ma x and r from Tableware presented in Figs. lll3l as functions of N part . Also 
there the lines of the best power approximations are depicted, 



X ~ Vp Qrt , X — Pmax i T, (23) 

with a scaling exponent k = 0.40 — 0.43 for p max and k = 0.23 — 0.33 for r. Note that the scaling is almost ideal, only 
slight deviations can be observed for the most central bins in some cases (p m ax and r for PHENIX at ^smn = 200 
GeV, Fig.^ and r for PHENIX at y / sjv r /v = 130 GeV, Fig.|3J). Another interesting point would be to give at least 

a rough notion of what the predicted volumes and V^ 2 ' mean. The simplest way is to compare them with the 
volume of the colliding nucleus, which is the gold nucleus here. Since the volume of the gold nucleus at rest (Va u ) 
equals about 1160 fm 3 (Va = 4/3 nr^A, tq = 1.12 fm), one can see that for the most central collisions at ^snn = 200 

GeV the volume at the beginning of the freeze-out, Jy* , is ~ 3V4„, whereas the volume when the freeze-out ceases at 

the central slice, Vf , equals w 5Va u - For the most central bin at ^s^n = 130 GeV this is w 1.9Va u and w 2.4Va u , 
respectively. As far as the lifespan of the system is considered the final moment of the freeze-out in the central slice, 

(2) 

t\ Q , has turned out to be about 10 fm for the most central bins. With the realistic assumption of the maximal 
possible value of the fluid rapidity, a™ ax , which will be explored in Sec. II V CI and for values of ajp aa: obtained there, 

(3) (2) 

the moment when the freeze-out ceases entirely, t\ , will be of the order of 10 t\ ' , that is of the order of 100 fm for 
the most central collisions. 

The extension of the model predictions for the invariant distributions to the low-p^ region (0.03-0.05 GeV for pions, 
0.09-0.13 GeV for kaons and 0.14-0.21 GeV for protons and antiprotons) agree well with the available PHOBOS data 
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TABLE I: Values of the geometric parameters of the model for various centrality bins fitted with the use of the RHIC final 
data for the pr spectra of identified charged hadrons fsl IsL l23 j I24I |25 | . Values of the thermal parameters are taken from the 
quoted references. 



Collision case 


Centrality [%] 


Npart 


Pmax [fm] 


r [fm] 


nmax 
P± 


V?>. [fm 3 ] 


tf.l H 


Vfi [fm 3 ] 


X 2 / ND F 


PHENIX at 


0-5 


347.7 


6.50±0.27 


8.23±0.23 


0.62 


2184.1 


10.5 


2782.2 


0.52 


Vsivjv = 130 GeV: 


5-15 


271.3 


5.99±0.21 


7.29±0.18 


0.63 


1641.1 


9.4 


2124.2 


0.46 


T = 165 MeV 


15-30 


180.2 


5.08±0.18 


6.34±0.15 


0.63 


1028.1 


8.1 


1317.6 


0.49 


= 41 MeV 


30-60 


78.5 


3.59±0.15 


4.81±0.13 


0.60 


388.9 


6.0 


485.4 


0.74 




60-92 


14.3 


1.68±0.19 


3.14±0.22 


0.47 


55.6 


3.6 


63.1 


1.32 


PHENIX at 


0-5 


351.4 


8.46±0.10 


8.84±0.08 


0.69 


3973.4 


12.2 


5497.0 


0.80 


Vsivjv = 200 GeV: 


5-10 


299.0 


7.99±0.10 


8.23±0.08 


0.70 


3302.6 


11.5 


4602.9 


0.61 


T = 155.2 MeV 


10-15 


253.9 


7.54±0.10 


7.67±0.08 


0.70 


2736.2 


10.8 


3837.1 


0.48 


fi B = 26.4 MeV 


15-20 


215.3 


7.11±0.10 


7.17±0.07 


0.70 


2275.5 


10.1 


3203.0 


0.48 


f20] 


20-30 


166.6 


6.45±0.09 


6.47±0.07 


0.71 


1689.5 


9.1 


2386.3 


0.58 




30-40 


114.2 


5.57±0.08 


5.63±0.06 


0.70 


1097.2 


7.9 


1544.0 


0.77 




40-50 


74.4 


4.68±0.07 


4.85±0.06 


0.69 


669.0 


6.7 


929.8 


1.05 




50-60 


45.5 


3.83±0.07 


4.16±0.05 


0.68 


383.9 


5.7 


521.6 


1.13 




60-70 


25.7 


2.99±0.06 


3.47±0.05 


0.65 


194.3 


4.6 


256.5 


1.41 




70-80 


13.4 


2.22±0.06 


2.78±0.05 


0.62 


86.3 


3.6 


110.4 


1.55 




80-92 


6.3 


1.71±0.06 


2.40±0.05 


0.58 


44.2 


2.9 


54.2 


1.40 


STAR at 


0-5 


352 


9.22±0.18 


7.13±0.06 


0.79 


3805.9 


11.7 


6222.1 


0.29 


Vsivjv = 200 GeV: 


5-10 


299 


8.40±0.17 


6.83±0.06 


0.78 


3029.0 


10.8 


4800.7 


0.27 


T = 160.0 MeV 


10-20 


234 


7.57±0.15 


6.33±0.06 


0.77 


2277.0 


9.9 


3548.3 


0.23 


fi B = 24.0 MeV 


20-30 


166 


6.50±0.14 


5.86±0.06 


0.74 


1557.9 


8.8 


2326.6 


0.30 




30-40 


115 


5.52±0.12 


5.37±0.06 


0.72 


1028.9 


7.7 


1476.2 


0.27 




40-50 


76 


4.66±0.11 


4.91±0.06 


0.69 


669.8 


6.8 


922.7 


0.27 




50-60 


47 


3.87±0.10 


4.40±0.06 


0.66 


413.3 


5.9 


550.2 


0.35 




60-70 


27 


3.07±0.09 


3.94±0.06 


0.61 


232.9 


5.0 


295.3 


0.46 




70-80 


14 


2.37±0.08 


3.32±0.06 


0.58 


116.8 


4.1 


143.4 


0.87 


BRAHMS at 


0-5 


357 


8.75±0.16 


8.38±0.11 


0.72 


4029.6 


12.1 


5825.3 


0.50 


y/s NN = 200 GeV: 


0-10 


328 


8.50±0.15 


8.08±0.10 


0.72 


3670.6 


11.7 


5329.1 


0.52 


T = 155.2 MeV 


10-20 


239 


7.52±0.13 


7.28±0.09 


0.72 


2584.1 


10.5 


3714.1 


0.46 


fi B = 26.4 MeV 


20-40 


140 


6.29±0.12 


6.20±0.09 


0.71 


1541.5 


8.8 


2194.7 


0.36 


[201 


40-60 


62 


4.42±0.12 


4.95±0.10 


0.67 


608.2 


6.6 


816.0 


0.61 



|28| . even protons and antiprotons agree within errors as one can see in Fig. 0] (but the m 30% overestimation of 
protons and antiprotons has been obtained for the lowest pr point of the PHENIX data, i.e. for pr — 0.65 GeV). 
The similar result was obtained also for the preliminary data but protons and antiprotons were much more overweight 
P|. In Fig.Q]the original PHENIX data Q have been modified so as to be more useful for the comparison with the 
PHOBOS data j^j. Namely, the PHOBOS measurements were done for the 15% most central collisions whereas the 
PHENIX ones for the — 5%, 5 — 10% and 10 — 15% centrality bins. Since the treatment of counts includes the 
averaging over the number of events in a given centrality bin and for the same run the number of events in the 15% 
most central bin should be equal to the sum of numbers of events in the — 5%, 5 — 10% and 10 — 15% centrality bins, 
the rough approximation of the hypothetic measurement done in the — 15% centrality bin would be the average of 
the measurements done in the — 5%, 5 — 10% and 10 — 15% centrality bins. Therefore, in Fig.0] the data points 
representing the PHENIX yields are such averages of the original values from Q . The model predictions have been 
averaged in the same way. 
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PHENIX 800 GeV 




a p max =0.80477N par t 

0.33107 



T =1.21059N 
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N 
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FIG. 1: Values of the geometric parameters of the model from the fourth and fifth column of Table [I] for PHENIX at 
^/sjvjv = 200 GeV. The lines are the best power approximations. 
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STAR 200 GleV 




a p max =0.77203N 
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t=1.80554N 
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100 200 300 400 



N 
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FIG. 2: Values of the geometric parameters of the model from the fourth and fifth column of Table|I]for STAR at ^/snn = 200 
GeV. The lines are the best power approximations. 



B. Estimations of transverse energy and charged particle multiplicity densities 



The results of numerical estimations of dN c h / drj\ m id divided by the number of participant pairs for various centrality 
classes are presented in Figs. 03 and for RHIC at ^Jsnn — 130 and 200 GeV, respectively. Additionally to the 
straightforward PHENIX measurements of the charged particle multiplicity density, also the data from the summing 
up of the integrated charged hadron yields 0, |2j| are depicted in these figures (more precisely, since the integrated 
charged hadron yields are given as rapidity densities, the transformation to pseudo-rapidity should be done, which 
means the division by a factor 1.2 here, see [27j]). Also in [g| the feeding of p(p) from A(A) decays is excluded. 
To diminish this effect, integrated p and p yields delivered in Q were corrected to include back the feeding. The 
correction was done by the division by a factor 0.65, which is the rough average of a p^-dependent multiplier used 
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BRAHMS 200 GeV 
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FIG. 3: Values of the geometric parameters of the model from the fourth and fifth column of Table [I] for BRAHMS at 
-y/sjvjv = 200 GeV and PHENIX at ^/snn = 130 GeV. The lines are the best power approximations. 
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FIG. 4: Invariant yields as a function of pt for RHIC at v / sjvjv = 200 GeV. The PHOBOS data are for the 15% most central 
collisions with the error bars expressed as the sum of the systematic and statistical uncertainties |28|. The corresponding 
PHENIX data are presented as the averages of the invariant yields for the — 5%, 5 — 10% and 10 — 15% centrality bins with 
no errors given. Lines are the appropriate predictions of the single-freeze-out model. 



by PHENIX Collaboration (see Fig. 8 and Eq.(5) in J8j). Generally, the agreement of the model predictions with the 
data is much better for RHIC at tJsnn = 200 GeV. For the case of ^/sJTn — 130 GeV, only the qualitative agreement 
has been reached. The main reason of the sw 15% underweight of dN c h/dr) seen in Figgis that the best fit estimates 
for the pt spectra of pions at the lowest px lies 29% on the average below the corresponding measured values. This 
concerns also the PHENIX case at ^/snn = 200 GeV but the underestimation of low-pr pions is smaller there (15% 
on the average). Since pions comprise about 80% of dN c h/drj and their low-p^ fraction can contribute even 30 % to 
the integrated yield 22], this explains the obtained underweight of dNch/dr] in spite of the very good quality of the fit 
for PHENIX at y^ijyjv = 130 GeV. The discrepancy between the directly measured dN c h/drj and dN c h/dr] expressed 
as the sum of the integrated charged hadron yields can be the next reason, especially for RHIC at ^/snn = 200 GeV 
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FIG. 5: dNch/dr] per pair of participants versus Npart for RHIC at ^/sjvjv = 130 GeV. The PHENIX original data are from 
[l^|. Also the recalculated PHENIX data from summing up the integrated charged hadron yields delivered in |2^l are depicted. 
The line connects the results and is to guide the eye. 

(see Fig.® this effect was notified in backup slides of |||). The discrepancy starts at mid-centrality and rises with the 
centrality. The another reason for the quantitative disagreement is that transverse momentum spectra are measured 
in limited ranges, so very important low-py regions are not covered by the data. To obtain integrated yields some 
extrapolations below and above the measured ranges are used. In fact these extrapolations are only analytical fits, but 
contributions from regions covered by them account for about 25 — 40% of the integrated yields |22( . It might turned 
out that these extrapolations differ from the thermal distributions supplemented by the distributions of products of 
decays. However, the extension of the model predictions for the invariant distributions to the low-px region agree 
well with the available data, cf. Fig.0] 

The values of cLEt / dr]\ m id divided by the number of participant pairs are shown in Figs.0and[S]for ^snn = 130 
and 200 GeV, respectively. The quality of the model predictions for dEx/di] is much better then for dN c h/dr). Again 
the worse agreement is for the case of ^/snn = 130 GeV. The STAR measurements need separate comments. The 
STAR data were taken not at mid-rapidity (77 = 0) as in the PHENIX case but at 77 = 0.5 [ly. But the pr spectra 
used for fits of the geometric parameters were measured at mid-rapidity, also in the STAR case |9j]. Therefore one 
should expect some overestimation of the predictions in comparison with the data on the whole. To remove this 
effect the original STAR data ^(| have been divided by a factor sin (0| ,7=0.5) ~ 0.887. As it can be seen from Fig. [SI 
the predictions and data agree with each other within errors besides two most central points. But even there, the 
discrepancy is w 17%. It should be also noticed that bigger values of transverse energy estimates for the STAR case 
are the consequence of slightly higher values of pt distributions measured by STAR Collaboration with respect to 
PHENIX measurements (it can be seen directly from careful comparison of spectra given in and ||). 

Values of the ratio (dEr / drj) / {dN c h/ drj) as a function of N part are presented in Figs. l9llTl As one can see, the 
position of model predictions is very regular and exactly resembles the configuration of the data in each case, the 
estimates are only shifted up about » 10% as a whole. For the PHENIX case, Figs. 1911 01 the shift is about 8% at 
\J snn = 130 GeV and about 14% at ^/snn = 200 GeV. This overestimation of the ratio can be explained, at least for 
more central collisions, by the observed discrepancy between the directly measured dN c h/drj and dN c h/dr] expressed 
as the sum of the integrated charged hadron yields. If the original data points are replaced by the recalculated data, 
such that the denominators are sums of the integrated charged hadron yields, then much better agreement can be 
reached for more central collisions. In the STAR case, Fig. ^2 as the experimental data at 77 = two modified versions 
of the original data 0] at 77 = 0.5 are depicted. The first version was done by the division of the original data by a 
factor sin (^^=0.5) ~ 0.887. In the second one, besides the above-mentioned rescaling, the values of dN c h/dr] measured 
directly at 77 = were put into the denominator. Up to the mid-centrality the model estimates agree pretty well 
with the second version of the data, but for the most central bin the overestimation of the ratio of about 13% has 
been obtained. 
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FIG. 6: dNch/dr) per pair of participants versus N par t for RHIC at ^/snn = 200 GeV. The original PHENIX data are from 
whereas the recalculated PHENIX data are from summing up the integrated charged hadron yields delivered in |§. Also 
the STAR data are depicted with no errors given as in the source paper |j| . The lines connect the results and are to guide the 
eye. 
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FIG. 7: dEr/dri per pair of participants versus N pa rt for RHIC at ^/JWn = 130 GeV. The PHENIX data are from [Eg . The 
line connects the results and is to guide the eye. 



C. Total charged-particle multiplicity 



Generally, for any expansion satisfying the condition da^ ~ on a freeze-out hypersurface, the total multiplicity 
of particle species i can be derived in the form (for the more formal proof see [loj') 



Ni = J d 2 p T dy d 2^ dy = J d2 PT dy J p^da^ / 4 (p • u) = J da J d 2 p T dy (p ■ u) f t (p ■ 
= J da J ^ (p ■ u) fi{p -u)= J da m(T,fi B ) = 7ij(T,/z B ) J da , 



(24) 
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FIG. 8: dEr/dr/ per pair of participants versus N part for RHIC at ^/snn = 200 GeV. The PHENIX data are from 0] but 
the original STAR data from llq ] have been rescaled to 77 = 0, see text for more explanations). The lines connect the results 
and are to guide the eye. 
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FIG. 9: (dE T /drj)/(dN ch /dri) versus iV part for RHIC at ^/IWn = 130 GeV. The original PHENIX data are from [T^. The 
recalculated PHENIX data are also depicted, here "recalculated" means that the sum of integrated charged hadron yields |22| 
have been substituted for the denominator in the ratio. The line connects the results and is to guide the eye. 



if the local thermal parameters are constant on this hypersurface. This means that under such condition the total 
multiplicity is a product of the density rii calculated for a static gas case and the hypersurface volume. Note that rii 
is not the primordial thermal density of particle species i but it collects also contributions from decays of resonances. 
Since the freeze-out hypersurface and the flow described in Sec. [n] fulfill the condition da^ ~ u M , the formula (|24l) is 
valid in this case as well. In practise the rapidity of the fluid element a\\ should not be unlimited but should have its 
maximal value a™ ax . Otherwise, as it can be seen further, the hypersurface volume and the total charged particle 
multiplicity would be infinite. Then, with the use of Eq. I|14|) one can express the hypersurface volume as 
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FIG. 10: {dE T /drj)/{dN ch /dTj) versus N part for RHIC at = 200 GeV. The original PHENIX data are from [l?]]. The 

recalculated PHENIX data are also depicted, here "recalculated" means that the sum of integrated charged hadron yields 
have been substituted for the denominator in the ratio. The line connects the results and is to guide the eye. 
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FIG. 11: {dE T /dri)/(dN ch /drj) versus N par t for RHIC at v /s^]v = 200 GeV. The original STAR data [3 have been rescaled 
to rj = by the division b y a factor 0.887. Also the recalculated STAR data are depicted, here "recalculated" means that the 
original values of dEr/dr/ |lfil ] have been divided by a factor 0.887 and values of dN c h/dr] measured at midrapidity [jj are taken 
for the denominator, however errors are not given in |9(. The line connects the results and is to guide the eye. 
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Finally, the total multiplicity of charged particles can be obtained: 
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One can see that the above formula has the form of the product of a volume and a density. This volume can be treated 
as an approximation of the size of the gas at the freeze-out (or an average volume of a hadron source) . Approximation 
- because it is the volume of a static fireball which would emit the same amount of charged hadrons as the freeze-out 
hypersurface described by Eq. Q with the superimposed limit in the range of an. Average - because for the static 
fireball the freeze-out happens at a given moment of time, whereas in the discussed model the freeze-out extends in 
time. This approximate volume will be called V s tat. from now on and equals 



V stat , = 2tt af ax rp 2 max . (27) 

For a™ ax the following reasonable assumption has been made: it is equal to the rapidity of leading baryons after 
the collision. In other words, it is assumed that the fluid which has been created in the CRR can not move faster in 
the longitudinal direction then fragments of a target or a projectile after the collision. Therefore it should depend on 
the centrality of the collision, since the more central the collision is, the higher degree of the stopping of the initial 
baryons ought to happen in principle. There are two limiting cases, the maximum stopping happens for the most 
central collision whereas if the centrality approaches 100% the stopping disappears. For values of a™ ax it means that 
for the most central collision a™ ax = (y) = y p — (Sy) and in the limiting case of the 100% centrality a™ ax = y p , 
where (y) is the mean net-baryon rapidity after the collision, y p the projectile rapidity and (Sy) the average rapidity 
loss [25j . And the last assumption is that a™ ax is a linear function of the centrality c, where c is a fractional number 
representing the middle of a given centrality bin, i.e. c = 0.025 for the — 5% centrality bin, c = 0.075 for the 5 — 10% 
centrality bin, etc.. Then ^^(c) has the following form: 



a T ax ^ = y^Mk- {1 - c) - (28) 

In the derivation of Eq. (|28|1 the most central bin has been taken as the — 5% centrality bin, since the estimate 
of the average rapidity loss (Sy) for RHIC at ^/snn = 200 GeV was done for this centrality class in |25j. The 
BRAHMS Collaboration reports (Sy) = 2.05 for RHIC at y/sJHf = 200 GeV (y p = 5.36) j^. Since for RHIC at 
Y / s]v]v = 130 GeV the average rapidity loss has not been given yet, it is assumed that (Sy) behaves linearly in y p 
between the SPS point (the 5% most central Pb-Pb collision at 158 GeV per nucleon: y p = 2.9 in the c.m.s. and 
(Sy) = 1.76 H|) and the BRAHMS one. This results in (Sy) = 2.0 for RHIC at = 130 GeV (y p = 4.93). The 

numerical estimates of n ch (T,/j, B ) for the static gas give n c h(T — 155.2 MeV, fi B — 26.4 MeV) — 0.36112 fm -3 and 
n ch {T = 160 MeV, \i B = 24 MeV) = 0.44839 fm~ 3 for RHIC at = 200 GeV and n ch (T = 165 MeV, [i B = 

41 MeV) = 0.5618 fm~ 3 for RHIC at ^snn — 130 GeV. Now the total multiplicity of charged particles N c h can 
be calculated with the use of Eqs. Ij26(l and (|28|l . The results presented as the total charged-particle multiplicity 
per participating pair versus N part are gathered in Fig.^| The predictions for STAR and PHENIX exhibit almost 
ideal centrality independence of the total charged-particle multiplicity per participating pair within the range of the 
PHOBOS measurement, i.e. N part ~ 60 — 360. For the BRAHMS case there is some deviation from this ideal 
behavior, but it is of the order of 10%. Additionally, the predicted values agree with the data within w 10%. 



D. Comparison with the statistical hadronization model 



Another very interesting point would be the comparison of the predictions of this model and the statistical hadroniza- 
tion model (SHM) of Ref. |2jj, from where values of the statistical parameters have been taken as an input to fit 
PHENIX and BRAHMS spectra at ^/sWn = 200 GeV in Sec. ITvTl An indirect comparison has been already done 
since good quality fits to the spectra have been obtained with the substitution of SHM statistical parameter esti- 
mates. This means that fits to the spectra with this model confirm values of statistical parameters obtained from fits 
to particle yields with SHM. The only left quantity which is independently obtained in these both models and could 
be directly compared is the volume of the hadron source. However not the global volume but the volume fraction 
dV/dy is predicted in |2jj. The volume fraction is defined by the relation |3Cj : 



dNi dV 

= ~T~ " n i ) (29) 

dy dy 

which in principle should hold for each measured particle species i. Note that dNi/dy depends on the centrality, so 
dV/dy does as well. According to Eq. (5) in |3fJ the volume of the hadron source is given by 



14 



d 



50 



40 



30 



20 - 



10 



PHOBOS 200 GeV 
PHOBOS 130 GeV 
UA5 (pp/pp) 300 GeV 
model: STAR 200 GeV 
model: BRAHMS 200 GeV 




model: PHENIX 200 GeV 
model: PHENIX 130 GeV 



100 



200 



300 



400 



" part 



FIG. 12: N ch per pair of participants versus N part for RHIC at ^s^v = 130 and 200 GeV. The PHOBOS data are from [IJ] 
and the pp/pp data point of the UA5 measurement is from Fig. 39. 5 in [l2T |. The lines connect the results and are to guide the 
eye. 



TABLE II: Comparison of predicted volumes of the hadron source for PHENIX at ^snn = 200 GeV. Values of dV/dy are the 
SHM estimates taken from Ref. For all volumes the unit is cubic Femtometers. 



Centrality [%] 


Npart 


Vstat. 


VsHM 


dV/dy 


0-5 


351.4 


13152 


8233-12349 


1920 


5-10 


299.0 


11279 


6899-10349 


1609 


10-15 


253.9 


9632 


5694-8542 


1328 


15-20 


215.3 


8250 


4961-7442 


1157 


20-30 


166.6 


6391 


3666-5499 


855 


30-40 


114.2 


4381 


2358-3538 


550 


40-50 


74.4 


2812 


1466-2200 


342 



Vshm = k — -2y p , (30) 
dy 

where k is a constant, k ~ 0.4 — 0.6 for RHIC and 2y p expresses the maximum rapidity range. Since this volume is 
in fact the volume of a static fireball, the most natural volume appearing in the single freeze-out model which Vshm 
can be compared with is V s tat. - the static equivalent of the freeze-out volume, defined by Eq. (|2Tjl . The detailed 
comparison of these both estimates of the volume of the hadron source for different centrality bins of PHENIX 
measurements at ^/snn — 200 GeV is given in Table ITU One can see that values of V s tat. and Vshm agree with each 
other qualitatively. Some overestimation of V s t a t. with respect to Vshm can be observed and this effect increases with 
the centrality. For the most central bin it is s» 6% relative to the upper limit of the Vshm range, whereas for the last 
listed one (40-50 %) the corresponding overestimation equals ~ 28%. 



V. CONCLUSIONS 



The single freeze-out model has been applied to estimate transverse energy and charged particle multiplicity densities 
as well as the total multiplicity of charged particles for different centrality bins of RHIC measurements at y / ijvlv = 130 
and 200 GeV. These three variables are independent observables, which means that they are measured independently 
of identified hadron spectroscopy. Since model fits were done to identified hadron data (particle yield ratios and pt 
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spectra) and the global variables are calculable in the single freeze-out model, it was natural to check whether their 
estimated values agree with the data. Generally the answer is yes, at least on the qualitative level. This means that 
the similar shapes of the functional dependence on the centrality of the model predictions have been obtained, only 
their normalizations differ w 10% from the experimental data. In particular, the most surprising result is the almost 
ideal independence of the total charged-particle multiplicity per participating pair on the centrality of the collision 
in the range of centralities covered by the PHOBOS data 11]. The range is N part ss 60 — 360, that is from the bin 
around 50% centrality up to the most central one. 

Even for dN c h/ 'dr\\ m id and dEx jdr\\ m id this analysis is nontrivial and does not mean simply a self-consistency check 
of various measurements. First, transverse momentum spectra are measured in limited ranges, so very important low- 
px regions are not covered by the data. For instance, for PHENIX the lowest pion point is at px = 0.25 GeV/c, the 
lowest kaon one at px = 0.45 GeV/c and the proton and antiproton data start at px = 0.65 GeV/c [Hl^. Therefore, 
to obtain integrated yields some extrapolations outside the measured ranges are used. In fact these extrapolations 
are only analytical fits without any physical reasoning, but, for instance, contributions from regions covered by them 
account for 30% of the yield for pions, 40% for kaons and 25% for protons and antiprotons for RHIC at ^snn = 130 
GeV 22]. On the other hand, a calorimeter acts very effectively for these species in the \ow-px range, namely pions 
with pt < 0.35 GeV/c, kaons with px < 0.64 GeV/c and protons and antiprotons with pt < 0.94 GeV/c deposit 
all their kinetic energy |Tsj . Second, it is impossible to check the consistency of the transverse energy data because 
not all stable hadron spectra are measured at midrapidity for each collision case. This mainly concerns neutrons and 
K\. Since the good predictions for the transverse energy density at midrapidity have been obtained (see Figs. 17180 . 
the present studies can be understood as an undirect proof that in these unmeasurable px regions spectra are also 
explicable by means of the thermal source with flow and decays. This has been also confirmed directly in the PHENIX 
case by the extension of the theoretical momentum distributions to the region of the very low transverse momentum 
covered by the PHOBOS detector (0.03-0.05 GeV for pions, 0.09-0.13 GeV for kaons and 0.14-0.21 GeV for protons 
and antiprotons) , see Fig. 0] 

Having compared the present results with the preliminary ones |l|, one can notice that substantial improvement 
has been achieved for the PHENIX case at -^saT/v = 200 GeV. In the previous analysis the behavior of the estimates 
of dEx I drj\ m id and dN c h / drj\ m id is very irregular in the region from peripheral to mid-central collisions and does not 
give the evidence for any clear dependence on the centrality. Here the dependence of predictions is exactly the same 
as the dependence of the data, only for dN c h / 'dr)\ m id the underestimation of at most 15% has been obtained. But the 
main success of this final analysis is the prediction of the scaling of the total charged-particle multiplicity with the 
number of participants for all considered cases, which failed in the studies based on the preliminary fits from 

As the next point, the possible description of the HBT radii within the single-freeze-out model will be briefly 
discussed (for a general review of the HBT interferometry see, e.g., |3l|). This problem needs some comments since 
the HBT analysis also gives information about the size of the hadron source. Studies of the HBT radii in the single- 
freeze-out model have been already done for the PHENIX data at ^/snn = 130 GeV 32]. It has turned out that the 
values of the estimated radii are about 30 % smaller then the measured ones. This discrepancy has been improved 
by the introduction of the excluded- volume corrections for the hadron gas, in the way as postulated in Ref. [33. The 
corrections give rise to the appearance of a new common scale factor, denoted as S~ 3 , in the invariant distribution, 
Eq. I|ll|) . This simply rescales both geometric parameters r and p max by the factor S, i.e. after fitting spectra t/S 
and p m ax IS have the same values as r and p m ax obtained for the point-like gas and listed in Table |U The value 
S = 1.3 has been found in Ref. what implies the increase of r and p m ax by 30 %. With these new rescaled 
parameters the predicted values of the HBT radii agree fairly well with the data [32]. Of course the detailed analysis 
of the HBT radii in the context of the single-freeze-out model for the final RHIC data at ^snn — 200 GeV should 
be performed and this will be the subject of further investigations. 

Finally, it should be notice that the statistical model of the present analysis corresponds to the chemical equilibrium 
case of the more general statistical hadronization model applied in Ref. [2(| (for a comprehensive review of the model 
see l^i) ). In some sense this analysis is complementary to those studies, because the estimates of statistical parameters 
from there have been taken as the input here and good quality fits have been obtained. Also the qualitative agreement 
between independent estimates of the volume of the hadron source done in both works has been obtained, c/. Table ITT1 
But as it was shown in Ref. [2fJ, the assumption of the strange quark non-equilibrium and further also the light 
quark non-equilibrium leads to more thorough description of particle yields for all accessible centralities of PHENIX 
measurements at ^snn = 200 GeV. Therefore such generalization of the present model should be done and tested 
with the data and this will be the subject of the next publication. 

To summarize, the single freeze-out version of a statistical (equilibrium) model fairly well explains the measured 
spectra of identified hadrons for various centralities and the centrality dependence of transverse energy and charged 
particle multiplicity pseudo-rapidity densities at mid-rapidity observed at RHIC. The fact which should be stressed 
again is that this model predicts the centrality independence of the total charged-particle multiplicity. Also predicted 
values of the total charged-particle multiplicity per participating pair agree with the measured ones within w 10%. 
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This is remarkable since geometric parameters have been fitted to spectra measured at midrapidity but the total 
charged-particle multiplicity deals with the whole rapidity range. 
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